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ABSTRACT 

We present first calculations of new long-distance contributions to B s — > 77 decay 
due to intermediate D s and D* meson states. The relevant 7 vertices are estimated 
using charge couplings and transition moment couplings. Within our uncertain- 
ties, we find that these long-distance contributions could be comparable to the 
known short-distance contributions. Since they have different Cabibbo-Kobayashi- 
Maskawa matrix-element factors, there may be an interesting possibility of observing 
CP violation in this decay. 
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A new era of experiments on rare B decays is about to dawn. The large data sets already 
obtained at LEP and CESR will be dwarfed by those provided by the e + e~ B factories, and 
by the hadron experiments HERA-B, CDF, DO and LHC-B. Among the rare B decays with 
particularly clean experimental signatures is B s — > 77, whose branching ratio presently has the 
experimental upper limit B{B S —> 77) < 1.48 x 10 -4 [p]]. Higher-sensitivity measurements of B s 
decays are not among the highest priorities of the e + e~ B factories, whose physics programmes 
are focussed initially on searches for CP violation in B& decays, but the hadron experiments 
cannot avoid being sensitive to B{B S — > 77) to levels several orders of magnitude below the 
present experimental upper limit [|IJ. 

The lowest-order short- distance contributions to the B s — > 77 decay arise from two sets of 
graphs: (i) box diagrams, and (ii) triangle diagrams with an external photon leg. These have 
been calculated 0-§], and yield, for m t « 175 GeV,[] B(B S 77) ~ 3.8 x 10" 7 , the next 
generation of experiments. Interestingly, the branching fraction can be enhanced substantially 
in extensions of the Standard Model such as a generic 2-Higgs scenario || . Emboldened by the 
fact that the short-distance QCD corrections are not too big 0-0], one may then even attempt 
to use such data as may be forthcoming as probes for new physics. 

However, it is not immediately obvious that this decay is dominated by such short- distance 
contributions. Estimates have been made of some long-distance contributions to the decay 
amplitudes, for example via intermediate charmonium states J / ip and ip' || . These were found 
to be very small, indicating that the short- distance contribution might be correct in order 
of magnitude. But, before concluding this to be the case, one must examine other possible 
long-distance contributions. 

In this paper we study such contributions due to intermediate D s and D* states via the 
diagrams shown in Figures |] and |^. These include loops of D s mesons alone, loops of D* mesons 
alone, and diagrams involving radiative D s — > D* transitions. The vertices are estimated using 
data on B — > DD, DD* and D*D* decays, the electromagnetic charge couplings of the D s 
and D*, and phenomeno logical estimates of the D* — > D s + 7 decay rate for the inelastic 
transition- matrix element. The imaginary part of the diagrams can be calculated easily while 
the determination of the real part requires the use of dispersion relations. 

We find that these long-distance contributions to B s — > 77 decay may be comparable to 
the short-distance contributions calculated previously, within the considerable uncertainties in- 
herent to the phenomenological inputs we use. We find no evidence that B(B S — > 77) decay 
should occur at a rate very close to the present experimental upper limit, but even this pos- 
sibility cannot be ruled out. Since the short- and long-distance contributions have different 
Cabibbo-Kobayashi-Maskawa matrix-element factors, the fact that they may be of compara- 
ble magnitudes offers the interesting possibility of observing CP violation in this decay mode. 
However, a detailed exploration of this possibility lies beyond the scope of this paper. 

At the quark level, the operators responsible for the B(B S — > 77) decay are b^^sF^F^ 
1 Over a wide range of the top-quark mass, the partial width grows linearly with m t [H. 
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and b^sF^Fw , where F is the electromagnetic fileld strength tensor, and F is its dual. The 
matrix element for the decay B s —>■ 7(^1^)7(92^) can then be parametrized as 
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where 1Z\ and IZ2 are the yet-to-be-determined hadronic matrix elements, and («(?) has been 
factored out in anticipation of the results to be derived later in the text. The partial width is 
then given by 
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It can be argued that the long-distance contributions should be dominated by two-meson in- 
termediate states. Now, the quark-level transitions with the least Cabibbo suppression are 
b — > ccs and b —* cud. The second process obviously leads to very dissimilar mesons which may 
rescatter into two photons only through higher-order terms. Consequently, we concentrate on 
the b — > ccs case, with the corresponding weak Hamiltonian being given by 



J, 
J" 



c7/x(! - 75)6 
S7 M (1 - 7s)c , 
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where V denotes the CKM matrix. The relevant intermediate states are then D+D~ , D* + D*~, 
D* + D~, etc.0. To calculate the B s 77 decay amplitude exactly, we need to determine the 
corresponding matrix elements of a non-trivial task. However, as a first approximation, 
we may assume naive factorization^ The if w k matrix elements are then expressible in terms of 

2 Since the contributions due to charmonia have already been investigated ||, we shall not consider these 
here. 

3 Based on experience elsewhere, we may hope that the associated errors are < 0(10%). 
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simpler quantities given by 
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where q = p — k in the last two terms, and m* = m£>*. The quantities V and Ai are related to 
the usual form factors through 
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Finally, we have 
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Within the factorization approximation,the relations in ([F]) can be interpreted as a parametriza- 
tion of the BD^D^ vertices. 

The form factors themselves are nonperturbative quantities. Apart from lattice calcula- 
tions 0, the best arenas for determining these are phenomenological approaches such as the 
BSW model fl(| or heavy-quark effective theory Some of these parameters are also mea- 
sured experimentally [Q, although a somewhat large spread persists in the data. Typically, 
though, we have /, /* ~ 200 MeV, f±, V ~ 0.6 and A\ )2 ~ 0.25. On the other hand, A3 is very 
small. 

As a next step, we need to determine the meson-photon couplings. Since we are dealing 
here with charged mesons, the simplest course is to assume minimal substitution. This then 
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fixes the charge coupling uniquely for each angular-momentum state of the (cs) system. Of 
course, the fact that these mesons are not fundamental particles means that we should ideally 
discuss a series of form factors, including charge radii as well as higher moments. As a first 
approximation, however, we neglect such aspects of the meson-photon coupling. 

The one-loop contribution due to the D s alone can then be expressed in terms of the 
diagrams in Fig. |l|. Computing these diagrams, the first one gives 




Figure 1: One-loop contribution to the B s — > 77 amplitude due to D s mesons alone. 
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In each of the above, the form-factors f± are to be evaluated at (p — k) 2 . 

We note that an exact calculation of the above integrals needs knowledge of the momentum 
dependence of the form factors. Furthermore, the integrals are formally divergent, and to 
calculate the real parts we would need to consider a cutoff scale that sets the limit of validity 
of a theory with 'fundamental' meson fields. Rather than attempt this, we take recourse to the 
optical theorem and calculate only the absorptive parts of the diagrams above. On using the 
Cutkosky rules, the sum of the three diagrams reduces to 
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Here = m 2 /m 2 B with m\ — m = niD s and 



m,£)*. The form of (R) is a testimonial 



to the fact that the three diagrams of Fig. [I] together form a gauge-invariant set. Thus, if D s 
were the only meson that could contribute, 7^1,2 would be determined by and 



ImKxiDs) = 0.39/ [f + {m 2 ) + 0.16 f-{m 2 ) 
lmJl 2 {D s ) = 
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The CP-violating form factor is thus identically zero. 



What about the real parts of the amplitudes TZi? These may be computed using dispersion 



relations ||14|| , and, for the case in hand, are seen to be smaller than the imaginary parts. 
Substituting for the form factors in the above and using (|3|), we then have 



Br(B s ^ 77) ~ 2.5 x 10" 



(10) 



This, by itself, is small compared to the short- distance contribution ||, though the interference 
term can be significant. 

What about the other long-distance effects? Of the ones calculated in the literature, B s — > 
07 — ¥ 77 has an amplitude somewhat smaller than the one calculated here, whilst B s — > 
J/ip<p —>■ 77 is seen to contribute only at the 1% level |§. We are thus in a situation where the 
2D S intermediate state may, in fact, give the largest long-distance contribution to the decay 
amplitude. It is interesting, at this stage, to compare the B s — > 77 case with that for K s — > 77 
decay. Whereas there it is the 2n intermediate state that dominates the decay amplitude |]I5| , 
the analogous contribution to the B s decay falls well short of the short-distance amplitude. 

We now turn our attention to the next higher state that can contribute to this process, 
namely the D*. As long as we neglect any D*D s '-f coupling, the additional diagrams involving 
the D* field are exact analogues of those in Fig. [I], and lead to 
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= 0.13 f*V(ml) 

Two new features confront us. One is that the CP-violating amplitude is non-zero, though 
small. More relevantly for our present purpose, we find lmTZi(D*) ~ 800 MeV, on the basis of 
which (|) then gives us 

Br(B s -A 77) ~ 6.5 x 10" 6 . (12) 

This, of course, is much larger than the short-distance contribution, and suggests that this 
particular decay should be observable at the very first run of the hadronic machines! 

A few objections could possibly be raised against this conclusion. For one thing, the validity 
of the factorization approximation as applied to decays into two vector mesons is not apparent. 
Thus, instead of using (||), it would be preferable to use the data to parametrize the B S D*D* 
vertex. Unfortunately, though, this exclusive mode is undetected so far. Once it is measured, 
it will be a straightforward task to reexpress our result as a prediction of the ratio of the two 
decay modes. 
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Secondly, we have, until now, neglected quite a few possible contributions to the process. For 
example, the contributions due to the higher excited states could in principle be comparable in 
magnitude, and might interfere destructively, thus reducing the total long-distance contribution. 
However, in the absence of extensive data on these states, it is almost impossible to estimate 
such contributions to any degree of reliability. We can only hope that this lack of information 
does not invalidate the results presented here. 

Finally, there is another class of contribution that we have neglected so far, involving the 
D*D S ^ vertex. Whilst this transition moment can, in principle, be calculated within a given 
model for the mesons, it is easier to work in terms of an effective Lagrangian, the relevant part 
of which can be expressed as 

. eV 
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We shall assume that T> is real, i.e. there is no absorptive part associated with this vertex. 
The new contributions are given by the diagrams of Fig. along with the crossed ones, and 
the resultant shifts STZi^ given by 
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Numerically, then, 

ImSKi = \V\ 2 {0.077/ [/ 
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One may, in principle, estimate T> from the partial decay width of D*: 
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T(D* S -> D sl ) = —\V\ Z mJ 1 - — 
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Unfortunately, this decay mode is not yet well measured. All we know is that r(-D*) < 1.9 MeV 
and that (D s + 7) and (D s + 71) are the only decay modes seen. This implies that 

\V\ 2 < 1300 Br{D* s -> D sl ) . 

A branching ratio of even 1% could then lead to an order of magnitude change in T(B S — > 77)! 
However, a stricter bound for T> can be obtained if one relates it to the corresponding form factor 



for the S = charm meson. Noting that T(D*) < 0.131MeV and Br(D* + -> D + ^) « 1.7% |3 



a relation analogous to (15) leads to 



\V d *d 7 \ 2 < 0.95 . 

A value of T> of this order obviously cannot negate the conclusions of (0). Examining 
closely, we see that the contribution due to a pair of on-shell D s exchanging a D* is actually 
much smaller than that of @. This result is similar again to the case of K$ —>■ 77 decay fl4| , 
where the analogous contributions due to the flavour-octet vector mesons are small. The other 
diagrams, where at least one D* s is on shell, can, however, compete with the 2D S contribution. 

To conclude, we have estimated the long-distance contributions to the B s — > 77 decay 
arising from charmed-meson intermediate states. We find that the 2D S contribution, by itself, 
is larger than the other long-distance contributions calculated hitherto in the literature ||. 
It is, however, still smaller than the short-distance amplitude. More interestingly, the 2D* 
contribution is much larger, and could enhance the branching fraction by more than an order 
of magnitude. This would lead to a very striking signal at the hadronic B factories. Moreover, 
if it is indeed comparable to the short-distance amplitude, the different Cabibbo-Kobayashi- 
Maskawa structures might offer interesting prospects for observing CP violation. 

DC acknowledges useful discussions with Ahmed Ali, Tariq Aziz, Leo Stodolsky and York- 
Peng Yao. 
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